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04 ' We investigate sohton colhsions in the Manakov model, which is a system of coupled 

nonlinear Schrodinger equations that is integrable via the inverse scattering method. Com- 
^ I puting the asymptotic forms of the general A^-soliton solution in the limits t Too, we 

pr . elucidate a mechanism that factorizes an 7V-soliton collision into a nonlinear superposition 

of (^) pair collisions with arbitrary order. This removes the misunderstanding that multi- 
OA ' particle effects exist in the Manakov model and provides a new "set-theoretical" solution 
to the quantum Yang-Baxter equation. As a by-product, we also obtain a new nontrivial 
relation among determinants and extended determinants. 

C I §1. Introduction 



In recent years, there has been a surge of interest in some systems of coupled nonhnear 
Schrodinger (coupled NLS) equations because of their relevance in nonlinear optics.—''—"^' Among 
^ ' such systems, this paper focuses on the following system of coupled NLS equations: 

a^ : 

\Q- i9t + + 2||q|| g = 0, g = (gi, 92, • • • , 9m)- (1-1) 

O ; 

■ Here ||qp = Q = SJLi kil^i where the superscript f denotes Hermitian conjugation. The 
, subscripts t and x denote the partial differentiation with respect to these variables. It is well 
' known that is a completely integrable system.'^^'^ We call the Manakov model, since 

the two-component (m = 2) case of was solved for the first time by Manako^^^ using the 
inverse scattering method (ISM). The extension of the ISM to the general m-component case is 
' straightforward.^' Nevertheless, the value of m is extremely important when we consider soliton 

■ solutions. The m = 2 case is merely a special case of the general m case. The most interesting is the 
^ , case in which the total number of solitons, say N , is equal to the number of components, m. Indeed, 

' in the N = m case, the coefficient vectors for the hyperbolic-secant-type envelope of solitons [see, 

■ e.g., ui in H2-12() ]. which we refer to after normalization as polarization vectors, are just sufficient 
to span the vector space C™", in which q exists. Soliton solutions in the other cases {N > m or 
N < m) are obtained from those in this case through a special choice of soliton parameters or the 
operation of a unitary transformation. In this paper, we consider the most general case, in which 
N and m are arbitrary positive integers. Then, the most interesting case, N = m,\s automatically 
included. 

Although the integrability of the Manakov model (jl-lj) has been established formally through 
application of the ISM, multi-soliton dynamics in the model remain to be clarified. There are two 

*^ E-mail: tsuchida@poisson.ms.u-tokyo.ac.jp 
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reasons for this. One reason is that the vector nature of (|1-H) . which supports the internal degrees 
of freedom of sohtons, leads to complicated behavior. Indeed, even a two-soliton collision is highly 
nontrivialj^'^^''^''^ as not only does a displacement of the soliton centers depend on the initial 
polarization vectors but also the polarization vectors themselves are changed. Therefore, in this 
model, the effect of an A^-soliton collision can never be written as the algebraic sum of the effects 
of pair collisions (at least, if we employ only the initial soliton parameters). This is quite different 
from the NLS case [i.e. with m = 1], in which the effect of an A^-soliton collision in fact can be 
written as the algebraic sum of the effects of pair collisions (in which the order of the pair collisions 
does not matter). The other reason is that Manakov gave a rather misleading description 
of an A^-soliton collision in Ref. |^. We quote the corresponding part, the first two sentences of the 
last paragraph of §2 in Ref. "l): 

Comparison of relations (17) and (18) indicates that an A'"-soliton collision does not, 
in general, reduce to a pair collision. This is clear, for example, from the fact that the ^j^^-, 
expression for 5^ contains with j > k, which depend on the initial parameters 
of all the remaining solitons. 
The relations (17) referred to here are given by 

■^^V = ^ n «ll^(CAf>Cn) } oP^{Cn,Ci,Si) . . .a^{CN,CN-l,S]^_-^)Sjj, 
{n<N ) 

st = \ll «ii(o, a) [ I n «n'(c- Cn) } «*(c, c.+i, • • • «*(C, Cn, s+) 

Kk>i ) Kn<i ) 

■a^iCi,Ci,S^)... a^(Ci, Ci-i, Sr^,)Sr, i = 1, 2, • • • , iV - 1, 
while the relations (18) are given by 

St = "iiHC2,Ci)a'^(C2,Ci,'S'5")5^, 
5+ = an(Ci,C2)«*(Cr,C2,.S+)5r. 

We now briefly explain the situation considered and the notation used in Ref.|4j|). The equations 

(17) of Ref. 1^ represent the solution for the collision of solitons (to which we refer as solitons- 1, 
2, • • • , A^), while the equations (18) represent that for two solitons. We note that the equations 

(18) are obtained from (17) by setting N = 2. The equations (17) were derived through analysis 
based on and intuition gained from the ISM. It is assumed that a soliton designated by a larger 
number moves faster along the x-axis. That is, soliton-i overtakes*^ solitons-1, 2, ■ ■ ■ , i — 1 and is 
overtaken by solitons-i + 1, i + 2, • • • , as time t passes from — oo to +oo. The velocity of soliton-i 
and its amplitude are determined by the complex parameter Q, which is time independent. The 
quantities Si are column vectors with two complex components, corresponding to the choice of 
m = 2 in Ref. The vector norm \\Si\\ [= (S] ■ Si) 2'^ determines the center position of soliton-i, 
while the normalized vector Sj/ISjl gives its polarization vector after the operation of Hermitian 

*' Throughout this paper, we use the term "overtake" if only the relative velocity is positive. Thus it can be used 
for head-on collisions, etc. 
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conjugation. The superscripts + and — denote the final state {t +00) and the initial state 
(t — > —00), respectively, an is a scalar function and d is a 2 x 2 matrix function. We do not give 
their explicit forms, which are not important in the following discussion. The superscripts T and 
* denote the operations of transposition and complex conjugation, respectively. 

Although the assertion (X) is somewhat ambiguous, the most natural and reasonable interpre- 
tation seems to be the following: 

Let us try to explain (17) by assuming that the solitons collide pairwise in ac- 
cordance with (18). Then, the first equation in (17) can be understood as follows. 
Soliton-A first overtakes soliton-A — 1 in its initial state, i.e. soliton-A — 1, which 
has not collided with other solitons. Next, soliton-A overtakes soliton-A — 2, which 
has not collided with other solitons. . . . Finally, it overtakes soliton-1, which has not 
collided with other solitons. 

Similarly, the second equation in (17) can be understood as follows. Soliton-i 
overtakes solitons-i — 1, i — 2, 1 in this order, none of which has collided with 
other solitons. Next, soliton-i is overtaken by solitons- A, A^ — 1, •••,i-|-lin their 
final states, i.e. those which will not collide with other solitons. 

If we attempt to diagram these events, we immediately encounter a contradiction. 
This indicates that an A-soliton collision cannot be described as a sequence of pair 
collisions, since the matrices a for different sets of arguments do not commute in 
general. 

The logic of the above interpretation is not mathematically rigorous, but it seems to be correct 
if the complex structure of (17) is taken into account. It might also be possible to interpret the 
assertion (Jft) in a different way. In any case, it appears that an A-soliton collision in the Manakov 
model ()1-1() does not reduce to a pair collision, and thus that some multi-particle effects exist in 
the Manakov model. However, in fact this is not true. 

The main goal of this paper is to clear up this misunderstanding. We explicitly demonstrate 
a mechanism that factorizes an A-soliton collision in the Manakov model ()1-1() into a nonlinear 
superposition of pair collisions. Here, we have used the term "nonlinear" to express the fact 
that the considered superposition is no longer additive. For the sake of definiteness, we explain 
in advance what we prove in terms of the Manakov notation, which also gives the definition of 
factorization in this paper. We first interpret the equations (18) as forming a nonlinear mapping 
with two complex parameters, /(C2;Ci)) which maps the initial state {S^2^,S^|f} into the final state 
{S^^, S^]''}. Then, we can use the mapping f{Cj,Ck) to evaluate in an A-soliton collision the effect 
of the two-soliton collision whereby soliton-j in a state Sj overtakes soliton-A: in a state Sk- For a 
given order of (^) pair collisions, we consider the corresponding composition of the (^) mappings: 
/(CjiCfe); N > j > k > 1. Then, regardless of the order of the pair collisions,*^ the composed 
mapping maps the initial state {S]^, S]^_^, • • • , S^} exactly to the final state {-S^^, S~^_-^^, • • • , Sf} 
given by the equations (17). 

To prove this factorization, we do not employ the Manakov results, equations (17) and (18), for 



*' Here, we are referring to the order in which neighboring sohtons collide pairwise. That is, unlike the m — 1 
case (scalar NLS), we do not consider virtual collisions between non- neighboring solitons. 
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the following two reasons. One reason is that, although it is ingenious and seems to be correct, the 
derivation of (17) in Ref. |1J) is neither very rigorous nor understandable to the reader not familiar 
with the ISM. The other reason is that the equations (17) are not tractable for our purpose. In this 
paper, we employ a more straightforward approach to obtain another formula for the asymptotic 
behavior of N solitons. We start from an explicit formula for the A^-soliton solution of the matrix 
NLS equation derived using the ISM in Ref . [T3|l . Through a simple reduction, we obtain an explicit 
formula for the general A^-soliton solution of the Manakov model (|1-1|) .^ To make the paper 
self-contained, we first set N = 2 and compute the asymptotic forms of the two-soliton solution 
in the t — > ^oo limits in our notation. These solutions define the collision laws of two solitons in 
the Manakov model, which are essentially the same as those given by the equations (18). Next, 
we consider the general N case and compute the asymptotic forms of the A^-soliton solution in the 
t — > ^oo limits. To express the polarization vectors appearing in the asymptotic forms concisely, 
we extend the definition of a determinant in such a way that the last column of an extended 
determinant consists of vectors. This determinant represents a vector defined in terms of the 
Laplace expansion with respect to the last column. We find a beautiful relation which casts the 
Hermitian product between such extended determinants into the form of a product of conventional 
determinants. Using this relation and the Jacobi formula for determinants, we prove that an A^- 
soliton collision in the Manakov model ()1-1() can be factorized into a nonlinear superposition of 
('^) pair collisions in arbitrary order. This reveals*^ the following two properties of the Manakov 
model: 

(a) An A^-soliton collision is composed of a nonlinear superposition of (^) pair collisions of certain 
order. 

(b) The composition of (^) mappings corresponding to pair collisions in an A^-soliton collision 
yields the same mapping for every possible order of composition. 

We remark that solitons are not mass points, and do not have compact support. Rather, they 
are structures with infinitely long tails. Even in a two-soliton collision, it takes an infinite time 
for the solitons to completely recover their own shapes. Taking this into consideration, it is more 
meaningful to understand this factorization conceptually than phenomenologically. It is unlikely 
that properties (a) and (b) are related directly. We note that proving (a) is equally difficult for 
every order of pair collisions and that proving (b) for arbitrary N reduces to proving it for = 3, 
that is, the following: 

(b') The composition of three mappings corresponding to pair collisions in a three-soliton collision 
does not depend on the (conceptional) order of pair collisions. 

The property (b') is called the Yang-Baxter property. The validity of this property means that the 
collision laws of two solitons in the Manakov model (|l-lj) give a "set-theoretical" solution to the 
quantum Yang-Baxter equation.^ To the best of the author's knowledge, this solution is new. We 
mention that another "set-theoretical" solution to the quantum Yang-Baxter equation is studied 
by Veselo^E^^ through investigation of the matrix KdV equation. 

A few interesting ideal^'^'^ have been proposed to explain in a general manner the pairwise 

*' To deduce property (b) from the factorization rigorously in the present setting, some discussion is needed. 
This is given in §5. 
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nature of soliton collisions in integrable systems. However, those ideas seem to be too intuitive in 
their present forms to complete a mathematically rigorous proof. For instance, it is not obvious 
that a pair collision is unaffected by the other solitons given only that they are far away, or that 
the final result is invariant when the order of multiple limits is changed. In addition, unlike this 
work, those works do not elucidate a nice mechanism of factorization. 

Finally, we would like to comment on the literature concerning the Manakov model (|1-H). Multi- 
soliton solutions of the Manakov model have already been obtained using the Hirota method^ ''^''^ 
(see also Refs. and I24[) for results obtained with another method). In this sense, although it 
is very useful, the explicit formula for the A^-soliton solution obtained in this paper may not be 
essentially new.*) The main contribution of this work is the elucidation of the pairwise nature of 
soliton collisions in the Manakov model. The results of this paper were obtained by the author in 
the summer of 2000 and presented at the autumn meeting of the Physical Society of Japan in that 
year. Very recently, he encountered some papers^^''^^''^'''**-' which pose the same problem (but 
do not solve it completely). In particular, in Ref. I27|). which actually appeared after the first sub- 
mission of this paper, Kanna and Lakshmanan gave a "proof" of the pairwise collision nature of a 
three-soliton collision. Unfortunately, the "proof" of Kanna-Lakshmanar^ is absolutely incorrect. 
Indeed, in addition to a few fatal mistakes, their "proof" as a whole is a typical example of circular 
reasoning. 

The remainder of this paper is organized as follows. In §2, we derive an explicit formula for 
the iV-soliton solution of the Manakov model.'^ In §3, we obtain the collision laws of two solitons. 
In §4, we compute the asymptotic forms of the A^-soliton solution in the limits t =Fco. In §5, we 
elucidate a mechanism that factorizes an A^-soliton collision into a nonlinear superposition of pair 
collisions and discuss the Yang-Baxter property. Section 6 is devoted to concluding remarks. 

§2. Explicit formula for the general AT-soliton solution 

In this section, considering a reduction of a formula given in Ref. 113(1 . we derive an explicit 
formula for the general A^-soliton solution of the Manakov model 

In Ref. 113(1 . under vanishing boundary conditions, we applied the ISM to nonlinear evolution 
equations associated with the generalized Zakharov-Shabat eigenvalue problem: 

" ^1 

Here, (" is the spectral parameter, I is the m x m unit matrix, and Q is an m x m matrix-valued 
potential function. The first two of the nonlinear evolution equations associated with ((2-l() are the 
matrix NLS equation, 

iQt + Qxx + 2QQ^Q = O, (2-2) 
and the matrix complex mKdV equation, 

Qt + Qxxx + 3{Q^Q^Q + QQ^Q^) = 0. (2-3) 

*' In any case, our formula has the advantage of compactness in its own right. 
**' Actually, the work of Park and ShirP^ is very similar to that of Steudel.l^ 



-Ki Q 
-Qt i^i 



^2 



(2-1) 
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We mention that integrable space-discretizations of (|2-2|) and (|2-3|) were found recentl}-^^' (see also 
the relevant work in Refs. i22j) and ,30,), 31), 32 j). The general A^-soliton solution of H2-2() or H2-3() is 
expressed aPJ 



Qix,t) = -2i(/I ••• I) S 

N 

where the mN x mN matrix S is given by 

N 



-1 



/ C7i(t)te-2iCi*- \ 
\ Cjv(t)"^e-2iC^^ / 



(2-4) 



Sjk = Sjkl - ^ 



tiici-cm-q) 



-Cj{t)^Ci{t), l<j,k<N. 



(2-5) 



Here, Q are discrete eigenvalues in the upper-half plane of C (ImCj > 0), each of which determines 
a bound state in the potential Q. The quantities Cj{t) are m x m nonzero matrices whose time 
dependences are given by 

C,(t) = C,(0)e^^^^, j = l,2,... ,iV, (2-6) 

for the matrix NLS equation ()2-2|) . and 

C,(i) = C,(0)e8i«^, i = l,2,... ,iV, 

for the matrix complex mKdV equation (|2-3|) . respectively. 

Let us consider a reduction of the A^-soliton solution of the matrix NLS equation to that of the 
Manakov model. We restrict the matrix Q to the form 



Q 



qi (12 









O 



(2-7) 



so that the matrix NLS equation H2-2|) is reduced to the Manakov model (jl-lj) . In this case, the 
matrices Cj(i)^ must have the same form as Q, from their definition,C3'^ 











(m) 



o 



i = l,2,-- - ,iV. 



(2- 



Conversely, if Cj(t)^ take the form (|2-8|) . Q{x,t) given by the formula (|2-4j) with (|2-5|) fits the form 
(|2-7|) . Then, the formula can be compressed into a compact form,^ 

N N 

q{x,t) = 2 j;(T-i),fce-2i?^c,(t), 
j=i k=i 
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where the N x N matrix T is given by 

^ e2i(G-c;)^ 

Thanks to 1)2-6^ and H2-8() . the time dependence of Cj(f) is given by 

c,{t)=e-^'^f'c,{0), j = l,2,--- ,N. 

The above set of equations gives a formula for the general A^-soliton solution of the Manakov model 
under vanishing boundary conditions. Let us rewrite this into a form convenient to investigate 
the asymptotic behavior. We first rewrite it as 

N N 

q{x,t) = 2^^(VF-i),fce-'[K^+«^)^+2(ci+Cf )*]c,(0), 
j=i k=i 

where the N x N matrix W is given by 
Next, we introduce the parametrization 

Cj = Cj + '^Vj eK, Vj >o), 

Cj{0) = 2r]je~"^Uj {aj e R, = 1), 

and employ the following abbreviations: 

Tj = -i[(0 - C;)x + 2(C| - Cf )t] = 2??j(x + 4Cjt), 
e, ^ (0 + C;)x + 2(C| + Cf)t = 2i,x + 4(^2 _ ^2)^^ 

In this way, we obtain the simplest formula for the general A^-soliton solution of the Manakov model 

N N 

q{x, t) = 2 ^ J^(C/-^),fce-'®'=Wfc, (2-9) 
j=i k=i 

where the N x N matrix U is given by 



Ujk = + E A,fcze-(^'+"')+i(®'-®^), 1 < j, A; < TV, (2-10) 

1=1 

with 

2r]i{uj ■ uj) 

(Q - Cfejly - Co j 



8 



If we set = 1 in the above formula, we obtain the one-sohton solution, 

q{x,t) = 277isech(ri +ai)e~'®i'Ui. (2-12) 

Therefore, we understand the significance of each parameter /coordinate as follows: 

2rjj : amplitude of soliton-j, 

— : velocity of soliton-j 's envelope, 

Tj : coordinate for observing soliton-j's envelope, 

Qj : coordinate for observing soliton-j's carrier waves, 

Uj : polarization vector of soliton-j (||itj|| = 1). 

To be precise, in the case of two or more solitons, the real polarization vectors are not invariant 
and change under soliton collision. The vector uj defines the bare polarization of soliton-j, which is 
realized when it becomes the rightmost soliton. This point is demonstrated below. In the following, 
we assume that all the soliton velocities are distinct, so that every soliton collides with all others. 



§3. Two-soliton collision 



In this section, we compute the asymptotic forms of the two-soliton solution in the limits 
t ^ =Foo, which define the collision laws of two solitons in the Manakov model 

We first write out the two-soliton solution given by (|2-9p with N = 2. According to (|2-lUp . the 
matrix U in this case takes the form 



U 



2i(i 



iize 



-{Ti+ai)+i{ei-ei) 



-(r,+a;)+i(©i-0l) 



1=1 



1=1 



1=1 



21/e 



-{Ti+ai)+i{ei-02) 



2r/2 



-(r;+Q!O+i(©i-02) 



Then, the two-soliton solution is given by 

2 



q{x,t) 



detC/ 



2r/2 



22ie 



-(ri+Qi)+i(0I-02) 



1=1 



E^ 

1=1 



21/6 



-{Ti+ai)+i{ei-e2) 



+ 



-,Ti+ai 



2r/i 



+ y^Aii;e" 



(Ti+ai)+i(e,-ei) 



1=1 



e-^®^«2 ) , (3-1) 



with 



det U = — h 



2r]i 2r/2 2r/i 



^ A22/e-(^'+"')+'(®'-®^) + 



.T2+a2 



1=1 



2m 



1=1 



_i_ g-(Ti+ai)g-{T2+«2) 

{/l,/2}={l,2} 



Aii/i X121 

A2I/2 ^22i 



(3-2) 
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Here, we have simplified the expression of det U using the relations 

0, / = 1,2, 



Alii Ai2z 

^211 A22Z 



(3-3) 



which can be proved straightforwardly. 
Next, we assume that 

6(= ReCi) >e2(=ReC2) 

and investigate the asymptotic behavior of q{x,t) as t ^ ^poo. This is accomplished by identifying 
the dominant terms in the numerator of l|3-l|) and its denominator (|3-2j) . We here note the relation 
n/rji = T2/t]2 + 8(6 - 6)*- 

In the limit t — > — 00, we have 

7"! To 
— < — . 

In this case, we have to consider separately the two regions (1^) and (2^) defined below. It is 
easily seen that q ~ in all other regions. 
(1~) finite ri, T2 ^ +00 

In this case, the dominant terms are those which contain the factor e^^. Then, using the 

relation Am = l/(2?7i), we obtain 



Ui 



2m 2)72 + 2r?2 

2r]i sech(ri + ai)e^^®^Mi. 



Aiiie-(^i+"i) 



(3-4) 



(2 ) Ti — 5- —00, finite T2 

Here, the dominant terms are those which contain the factor e"'^^. Then, we obtain 

2 [(A221 - A2ii)e-(-i+°i)e-'®2wi + (Am - Ai2i)e-(-i+"i)e-'®2n2] 



0^2 +"2 

2r;2 



Aiiie-(^i+"i) + e-(^i+°i)e-(^2+a2) ^^^^ 

={1,2} 



Aii/i 
A21Z2 



Al2i 

^221 



(3-5) 



In terms of i;^>i2 defined by 



-2(^12 



1 



A111A222 



E 



Aii/i 



Al2«i 



{/i,«2}={l,2} 

we can rewrite the asymptotic form (|3-5|1 as 

A221 



AoiZo Ao 



q ~ 2?72 sech(r2 + 02 + 4>i2)e x e' 



211 



A 



111 



\22«2 



Wl + 1 



(3-6) 



A 



121 



A 



111 



U2 



(3-7) 



Here, i;^i2 is always taken as real, since (|3-6|) can be rewritten as [cf. (|3-3|) and (|2TH1 

1 



-2cf>i 



A111A222 



Am + A112 A121 + A122 
A211 + A212 A221 + A222 
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(2r/i)(2r?2) 



.UiUl .Uiul 

.U2U\ .U2-ul 

1 — 1 — 

'Ci-Q C2-C2* 



2r)i 



Ci-Ci* Xi-C: 



2m 



C2-cr ^C2-C2* 

2' 



C1-C2 



Ci - C2* 



fl , (Cl-Cr)(C2-C2) 

1 ici-ai^ 



til • ''^2 



(> 0). 



In the limit t — > +00, we have 



??1 f?2 



In this case, we have to consider separately the two regions (2"'") and (l"*") defined below. It is 
easily seen that g ~ in all other regions. 
(2"'") Ti — > +00, finite T2 

In this case, the dominant terms are those which contain the factor e"^^ . Then, using the 

relation A222 = l/(2?72), we obtain 

2^e-e^.2 



2j7i 2»72 2r;i 

= 2r?2 sech(r2 + Q2)e^^®2'U2. (3-8) 

(1+) finite Ti, T2 — > —00 

In this case, the dominant terms are those which contain the factor e"'^^ . Then, with the help 
of H3-6() . we obtain 



q 



2 [(A222 — ^212)6" 



-(r2+a2)p-i© 



e-'^i^i + (A112 - Ai22)e- 



-(T2+a2)p-i©i 



U2 



2^r^A222e-(^2+a2) + e-(^i+"i)e-(^2+"2) i,} 

={1,2} 



A2I/2 -^22«2 



2?7i sech(ri + ai + (pu)^ x e' 



1 A212\ , / A112 — A122 , 

1 - T ] Ui+ I ^ U2 



A 



222 



A 



222 



(3-9) 



Taking the sum of (|3-4|) and (|3-7|) . or (|3-8|) and (|3-9|) . with a slight simplification, we arrive at the 
following theorem. 



Theorem 3.1. The asymptotic forms of the two-soliton solution of the Manakov model are 
as follows {see also Fig. 
as t ^ —00, 

q ~ 2r]i sech(ri + ai)e^'®^Wi + 2772 sech(T2 + 02 + (?!)i2)e~'®2it{i}^2; 

as t ^ +00, 

q ~ 2r/i sech(ri + ai + '/'i2)e~'®^ti{2},i + 2t/2 sech(T2 + a2)e~'®2u2- 
/:fere 012 a?^c? ti{i},2i '^{2},i o-re given by 



-2(t>i 



C1-C2 



Ci - CI 



1 + 



(Ci-cr)(C2-a) 

ICi-ClP 



Itl • 
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and 



u 



{1},2 



u 



{2},1 



Ci 



a 
cr 



C2-cr 



U2 



Ui 



Ci - cr 



Ci 

C2 



C2* 

CI 



C2 - cr 



U2-U{ \ Ui 



Theorem 13.11 defines the colHsion laws of two sohtons in the Manakov model, which we use in fJHl 
to factorize an A^-soliton collision into pair collisions. Here, we mention some important properties 
of the collision laws: 

• A two-soliton collision changes neither the amplitudes of the solitons nor the modulus of the 
Hermitian product of the polarization vectors. In fact, recalling that = ||tt2|| = 1, we can 
prove by direct computations that 



{1},2| 



U 



{2},l| 



1, 



Ul ■ u 



t 

{1},2 



U 



{2},1 



U 



Ul ■ U2 



u 



{2},1 • "{1},2 



Although we omit here the tiresome proof, the most important relation || Wj^} 2 1| = ||it{2},i || = 1 
is shown in ^ in a more general context. This relation shows that the collision is elastic if we 
observe it with conserved density, ||gp = Yl^=i kiP- 

As a result of the collision, the polarization vectors rotate nontrivially on the unit sphere in 
C™. Thus, if we observe the collision with respect to each component q^, it appears as if it 
were inelastic. 

We have expressed (pi2, '^{i},2 ""{2},! terms of ui and U2- Then, the collision laws are 
symmetric with respect to interchange of the subscripts 1 and 2. This form of the collision 
laws is very useful in studying the factorization of an A^-soliton collision into pair collisions. 
For any fixed unit vector Ui, we can invert the mapping U2 ^ ''^{i},2 using the following 
relation for the projection operator u\ui (cf. {u^ui)"^ = u[ui): 



Ci - C2* ' 7 V ' cr - C2 



I. 



Then, we can express 0i2, U2 and U{2},i in terms of ui and Wj^} 2 as Manakov did in Ref. 0)): 



-2</<12 



U2 



C1-C2 



Ci - C2* 

.-4>i2 Ci - C2 

^1 ^2 



1 



C2-C 



(Ci-cr)(C2-a) 
IC1-C212 



Ul ■ uj^j 2 



, Ci - cr f t, 

U{1},2 + 7737^ ( "{1},2 Ul 



-(Aio C2 ~ Cl , 

U{2},i = e "^12 7 <; Ul 



C2 - CI 
C2-C1 



Ul ■ u 



{1},2 "{1},2 



(3- 10a) 
(3-lOb) 
(3-lOc) 



Owing to the lack of symmetry under the exchange of 1 and 2 (even after a change of the 
vector notation), this new form of the collision laws is not useful in studying the factorization 
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problem. On the other hand, it is of prime importance in exphcitly showing that a two-soliton 
cohision can be expressed as a mapping from the initial state to the final state. Moreover, 
using H3-1U() . we can easily check that if ||mi|| = ||it{i},2|| = 1; then \\u2\\ = 1. This verifies that 
for any unit vector ui, the mapping U2 ^ ^{i},2 Theorem 13.11 is a bijection on the unit 
sphere. We use this fact in §5 to prove the Yang-Baxter property of the mapping (|3-1U() . 



q ~ 2r]i sech(ri + ai + 0i2)e '®^tt{2},i q — 2772 sech(r2 + 0^2)6 



soliton-1 



soliton-2 



A 



t 



/ \ 

; interaction j 
\ range ; 



soliton-2 



soliton-1 



q ~ 2r/2sech(r2 + a2 + (/>i2)e '®2tt{i},2 Q — 2r/i sech(ri + ai)e 



Fig. 1. Twosoliton collision. 
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§4. Asymptotic behavior of the AT-soliton solution 



In this section, we compute the asymptotic forms of the A^-sohton solution in the hmits t — > ^oo 
and simphfy them as much as possible. Extensions of some techniques applied to the KdV equation 
are used.^^''^^'''^^' 

We first rewrite the A^-soliton solution (|2-9j) before considering the limits t — > =Fc«. We use the 
tilde to denote cofactors. For instance, the cofactor Ukj is obtained by deleting the /c-th row and 
the j-th. column from the determinant of U and multiplying it by (—1)'^ . Using the definition of 
U given in H2-1U() and the multilinearity of determinants, we can rewrite (|2-9() as 



q{x,t) 
2 

detU 
2 

detU 
2 

detU 



N N 



j=i k=i 



[{Uii + ■■■ + UiN)e-'^'ui + ■■■ + {Uni + ■■■ + UNN)e 



ie. 



detU 



1 1 •• 

Uni Un2 

N-l 

E E 

n=0 2<ji< - <jr„<Af 



1 

U2N 

Unn 

N 

n 

. k=2 



wi H + 



Uu 



'UN 
U12 



Un-11 Un^12 

1 1 



2% 



N 



h,--,ln = l 



X Q-in„+ai„)+K^i„-&jn) 



1 



1 



A 



A. 



UiN 

Un-in 
1 



N-l 

+ - + E E 

n=0 l<ji< -<in<Af-l 



^1 ^Tk+a^ 



n 



. k=l 



2% 



1 



N 



^ g-(r,^+aii)+i(e,j-0,j) 



1 1 



e-'^'^UN 



(4-1) 



Similarly, we can rewrite the determinant of U as 



detU 
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N 

n 

k=l 



2% 



N 

ii=i 



N 

n 

1 k=l 



\ 



2% 



N 



i<ii<i2<A^ 



AT 

n 

1 fc=i 



2% 



h=i 



N 



h,l2 = l 



J2 



-iT,^+a,^)+i{ei^-ei) 



-(f"iiV+"i^r)+i{©i^r-0iv) 



Alii, 



Al7V/i 



N 



E E 

n=0 l<ji<-- <jn<N 



N 

n 

1 fc=i 



2% 



AT 

E 



X > e 



'(Tii+«ii)+i(e,^-e,,) . . .p-(r,„+Q,,j+i(e,„-e3„) 



■^iiji'i 



A 



•jnjlln 



Here, we note [cf. the definition of Xjki H2-ll() ] that the quantity 



Ci-cp iCi-cp)iCi-Q){Ci-q) 

is invariant under interchange of the subscripts j and j' . Thus, we have 



Cl - Cy 



^jkl 



ci-q 



■ Ai'/ 



■j'ki 



0. 



(4-2) 



This shows that if uji-uj 7^ or Uj-uj ^ 0, the two vectors [Xju, Xj2i, ■ ■ ■ , Xjni) and {Xj'u, Xji2i, ■ • • , Xj' 
are hnearly dependent. In the case that uj' ■ uj = Uj -uj = 0, according to (|2-11() . both vectors be- 
come zero. Therefore, the determinants in ()4-l|) or H4-2p contribute only if /i, • • • , are all distinct. 
This fact is a generalization of the relations (|3-3p . 
Next, we assume that 

6(=ReCi) >6(=ReC2) >--->Cn{= ReCTv), 

and investigate the asymptotic behavior of q{x,t) in the limits t ^00. This is accomplished by 
identifying the dominant terms in the numerator of (|4-1() and its denominator ()4-2[) . We here note 
the relations Tj/r]j = Tk/r]^ + 8(^j — ik)i- 



In the limit t 



-00, we have 



Ti T2 tm 
— < — <•••<—. 
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In this case, we have to consider the following N regions (l^)-(A^^) separately. It is easily seen 
that q ~ in all other regions. 
(1~) finite ti, T2, - ■ ■ ,tn ^ +00 

In this case, the dominant terms are those which contain the factor e'^^"' ^'^^ . Then, using 

the relation Am = 1/(2771), we obtain 

2 (nl2 ^) e-'^^^^i 



k=i -^^r + \nk=2^2^) e ^ 1+ ^^Aiii 

2r]i sech(Ti + ai)e^^®^iti. 



(4-3) 



(n ) Ti, • • • ,T„_i ^ -00, finite T„, t„+i, • ■ ■ , tat ^ +00, n = 2, •••,A^-1 

Here, the dominant terms are those which contain the factor e"'^^ Tn-i+rn+i-\ hrjv_ Then, 
those in the numerator of ()4-l() are 



n / N 

n 

j=l \k=n+l 



2% 



E 



-(n+ai) _ _ _ -(r„„i+a„„i) i(ej-e„) 



1 



={l,-,n-l} 
Alnii 



1 



Aj+lnij 



while those in the denominator (|4-2() are 



■J' 



n^) ^ 

={l,-,n-l} 



^(n^) E 



g-(Ti+cn) . . . g— (t„_i+q:„_i) 



-(ri+ai) _ _ _ „-(T„+a„) 



Aiizi 

An-lli„_i 
Allii • • • Ainjj 



Aln-lJi 
An— 1 71— 1 



A, 



Armi„ 



--{l,-,n} 

As a natural extension of ()3-6|) . we define (jji^i^.-.i^ for distinct positive integers 11,12,- " > V '-'y 



1 



A7-I 7,1 * * * Aj, 



E 



{h,-,lp} 

= {il,---,ip} 



A, 



A,; 



(4-4) 
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We should note that 4>iii2---ip is symmetric with respect to permutations of the subscripts 
h,'i'2i • • • ,'i'p- We prove below that (pi^i^-.-ip can always be taken as real. In terms of (pi^i^-^-ip, 
we can express the asymptotic form of q in this region as 



q ~ 2r]n sech(rn + + 



'12-ra 



!'i2-n-i)e"^®"e'^i2...n+<Ai2...n-i 



viii • 



={l,-,n-l} 



1 



1 

^nnl„-i 



It, 



(4-5) 



(A^~) Ti, • • • ,riv_i ^ — cxD, finite tat 

Here, the dominant terms are those which contain the factor e~'^^ "^iv-i^ "With calculations 
similar to those in the case {n~), we obtain the asymptotic form of q given by (|4-5j) with 
n = N. 



In the limit t — > +oo, we have 



Tl T2 Tat 

— » — »•••> — . 

Vl V2 VN 



In this case, we have to consider the following N regions (A^^)-(l"'") separately. It is easily seen 
that q ~ in all other regions. 
(A^+) n, • • • , TN-i +00, finite tn 

In this case, the dominant terms are those which contain the factor e'^^"' '"'^Jv-i. Then, using 

the relation Xnnn = l/(2??Ar), we obtain 



2r/fe \^llfc=l 2% 

2rj]sf sech(Tiv + aAr)e~^®'^ujv. 



(4-6) 



(n+) Tl, • • • ,T„_i ^ +00, finite r,i, t„+i, ■ ■ ■ , tat ^ -oo, n = 2, •••,Af-l 

In this case, the dominant terms are those which contain the factor e'^^"' ^'^n-i-r^+i- 
Then, those in the numerator of ()4-l() are 



-TN 



N /n-1 

n 

j=n \fc=l 



2% 



E 



g-(Tn + l+a7i + l) . . . g- 



(r]v+ajv)gi{ej-e„) 



{il,--- ,ijv-n} 

={n+l,--- ,Af} 



18 



A 



1 



A 



^j-lNlj^„ 
1 

'j+lNlj^„+i 



while those in the denominator (14-21) are 



■J' 



n 



E 



-(Tn+l+On + l) . . . p-(TJV+Q:jv) 



={«+!,■■■ ,N} 



I'n-l 



Kk=l / {h,-,lN-«+l} 

={n,- - ,N} 



nl 



N-n + l 



XnNh 
Aatat/jv-u+i 



In terms of 4'iii2--ip defined by (|4-4|) . we can express the asymptotic form of q in this region as 
q ~ 2rjn sech(r„ + a„ + </)„„+i...7V - </'n+i„+2...7v)e-'®"e^""+i-^+<^"+i"+2-'^ 

Ajinii ■ ■ ■ XnNh 



1 



N 



An+ln.+ln+l ' ' ' ^nnn 



E E 

={n+l,-,N} 



1 



1 



^NNU 



u,. (4-7) 



(1+) finite Ti, r2, • • • , tat ^ -oo 

In this case, the dominant terms are those which contain the factor e""^^ With calcu- 

lations similar to those in the case (n"*"), we obtain the asymptotic form of q given by (|4-7|) 
with n = 1. 

We now simplify the above asymptotic forms by using the fohowing definitions: 



Ck - q 



djk 



(4-8) 



Ck - q 

According to the definition of Xj^i H2-ll() . we have Xj^i = 2r]iCkidji. Then, we can rewrite the 
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definition of [ (|4-4|) witli p ^ n] in a factorized form [cf. the paragraph below (|4-2|) ]: 



1 = 1 



l\ — "^l)'" fin — ill'" i'^n 



2r}} Cj.i dj 1 •■• 2?7/ c,- / dj i 



1=1 



Hill 



HiX„ 



di. 



di 



(4-9) 



Here ii, i2, • " " j are distinct positive integers. For any nonzero vector (yj-^, ■ ■ ■ ,yi^), we have 



/ d 



{vh ) ■ ■ ■ ) y«n ) 



\ d 



V yl I 



E 



j,k=ii,--- ,in 
f oo 







E 



iU 



]=ll,--- ,in 



dz 



> 0. 



Thus, the eigenvalues of the underlined Hermitian matrix are all positive. This proves that the 
second term on the right-hand side of ()4-9p is positive. Considering the special case in which all 
the vectors , • " " ) are identical, we can prove the same for the third term in (|4-9|) . Therefore, 
the right-hand side of ()4-9() is positive, and 4'iii2---in can always be taken as real. In the same way 
as for (|4-9p . we can rewrite the second line of H4-5p or 1)4 -Tf) in the following factorized form: 



Ai.i 7,1 i-i ' * * A 



E E 

i=i {Zi, ■■■,/„-!} 

={ii,-- ,in-l} 



A 



iiiih 



X 



ilinh 



^ij^iiil 



1 



^ij-linl 

1 



J-1 



A 



A, 



'^j-\-l'^n^j 



Ui 



n-l 



11(2^.) xE E ••• 



E 



'^VhCiiiid; 



iih 



'^miCi^iAiii 



1 

'^%^hijdi,+ii, 



1 

2^ij Ci^i-di.^-^^i- 



'^'ni„-l'^ill„-ldi„l„_i ••• '^Vln^l'^inln-ldinln-1 
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1=1 j=i 





d. 



n-1 



-till 



1 



diiin-i 

1 



(i,- , 

*ntn — 1 



Ci 



1 



^in^n— 1 
1 



^^n^Ti — 1 
1 



diiii 



dinii 



di„i„-i '"■in 



(4-10) 



The last determinant, which contains vectors in its last column, represents a vector defined in terms 
of the Laplace expansion with respect to the last column. 

We can simplify the asymptotic forms further by noting some relations between the conventional 
determinants in (|4-9|) and ()4-10|) . We have the following lemma: 



Lemma 4.1. The following equalities involving determinants hold: 



-Jill 



"lull 



1 = 1 



^U?=liCin -Cil) 

nr=i(o„-a) 



ilii 



n-1 



A. _ A* n 



1 



^inin — 1 
1 

'^n-lil 
Cj„_l j„_l 



(4.11) 



(4-12) 



(4-13) 



Proof. We can prove (|4-11|) by subtracting on the left-hand side the last column from each of the 
other columns and using the relation 

/■^ 



Cjk — Cnk = A _ A* Cjfc. 

Similarly, H4-12p is proved by subtracting on the left-hand side the last row from each of the other 
rows and using the relation 

Cn - Ck 



Cjk C-jji 



Cn-C 



:Cjk- 
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4-13|) is a direct consequence of (j4-12|) and (|4-ll|l 



□ 



Taking the sum of and (jl^I) (n = 2, • • • , A^), or (jlH) and ^1^) {n =!,■■■ ,N -1), with 
the help of ()4-9() . H4-1U() and Lemma l4.ll we finally arrive at the following proposition. 

Proposition 4.2. The asymptotic forms of the N-soliton solution of the Manakov model (j i- j|l are 
as follows: 
as t — > — oo, 

TV 

q ~ ^ 2r/„ sech(r„ + a„ + ,„_i},„)e~'®"M{i,... ,„_i},„; 



as t —i- +00, 



n=l 



AT 



q-^^Vn Sech(r„ + a„ + (/){„+!,... ,Ar},n)e ^;v},n- 
n=l 

i/ere (/'{j^,... "^^jn,--- ,i„_i},i„ ore defined for distinct positive integers ii, - ■■ ,in-i,in by 



n-1 

n 

1=1 





diiii 


diiin 






di„ii 


dinin 






diiii 


diiin-i 


dinin 










di„-iii 


in — lin — l 



>0), (4-14) 



and 



u 



{il,--- ,in-l},in 



i = l {h,-,ln-l} 
= {il,--- ,in-l} 



^iiilii • • • ^i„-iin-iin-i 



\iiili 



^ij-iiil 



A 



A 



A 



^71^1 in — 1 



ij — 1 in^j — 1 

1 

^inin^n—1 



Ui 



Hnin-1 '^in 



diiin-1 



din-lil ' ' ' di^_-^i^_j^ 



(4-15) 
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When the set {ii,--- ,in-i} is empty, the definitions i4'Ml o,nd ( |^- j5| ) should read e = 1 
and it{ y^i = Ui. 

We prove in the next section that it{ii,... is always a unit vector, i.e. - || = 

1. This ensures that an A'^-sohton collision does not change the amplitudes of solitons. The 

vector ,n-i},n gives the polarization vector of soliton-n before an A^-soliton collision, while 
... „ gives that after the collision. Using the definition 

%i,-,in-i},i„ = 2??inSech(ri„ + ai„ +0|i^,...,i^_^},iJe-'®^"W{i,,...,i„_,|,i„, (4-16) 

we can diagram the asymptotic behavior of the A^-soliton solution in the simplest way (see Fig. I^J. 
We should note that Qji^,... is symmetric with respect to permutations of ii, • • • ,in-i- The 

last subscript i„ denotes the soliton's number, which is, of course, time independent. The signifi- 
cance of the other subscripts, ii, • • • ,?n,-i, in { } is clarified in the next section. 
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^{2,...,iV},l 



A 



t 



<?{n+l,...,7V}, 



n 



<7{},iv 




n 




interaction 
range 




f 



n 



^{l,...,Ar-l},iV ^{l,...,n-l},n 



Fig. 2. Asymptotic behavior of the A'^-soliton solution. 



24 



§5. Factorization of an AT-soliton collision into a superposition of pair collisions 



In this section, on the basis of the collision laws of two solitons presented in ^ we examine the 
asymptotic behavior of the A^-soliton solution obtained in ^ We conclude that an A^-soliton colli- 
sion in the Manakov model can be factorized into a nonlinear superposition of pair collisions 
in arbitrary order. 

We first prove a lemma needed later to compute the Hermitian product of ... j and 

■"{ii,--- ,i„-i},k- 



Lemma 5.1. For any set of unit vectors Ui^ , , Uj,Uk and du defined by ^4-^ , the following 
equality holds: 



diiii 



Ck - Q 



dji„_i Uj 
diiii ■ ■ ■ di^i^_^ 

di„_iii ■ ■ ■ di^_-^i^_-^ 



t 



di„_iii ■ ■ ■ di^_^i^_^ '"■in-l 



d 



d-jin-i djk 



(5-1) 



Proof. In the proof of this lemma, we use D to denote the last determinant in (jS-ip : 

diiii ■ ■ ■ C^jii„_i dij^k 



D = 

We express minor determinants obtained by deleting one row and one column of D as 



dji„_i djk 



D 



k 



diiii 



diiin-i 



d 



d 



Hi_iln-1 
H + li-n-l 



D 



3 
il 



diiii 
din-ih 



di„^iii 
djii 

diiii_i diii^^-^ 



di„_ii„_i 



d 



diii„_i diik 



di„_iii_i di^_iii_^i ■ ■ ■ di^_ii^_i di^ik 
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D 



J 
k 



di\in-i 



d 



di. 



Using these abbreviations and the Laplace expansion of determinants, we can rewrite the left-hand 
side of (|5-1|) as 



( n-1 






3 
k 










l.h.S. = <^ 


Ip 
k 






. h . 


ul +D 


3 
k 


[p=l 











u 



n—1 n~l 

p=l q=l 
n-1 

9=1 
n-1 

p=l 



D 



(c.,-c;) + (c;-c4) 



3 
k 



D 



D 



3 
k 



(a-c;) + (c; -c) 



3 
k 



n-1 

q=l 
n-1 

p=l 



—d 



'jk 



(Ci,-C) 



n-1 
p=l 



k 



3 
k 



k 



n-1 
q=l 



d^,^, + D 



3 

k 



dipk 



+ D 



n-1 



3 

k 



iCk - q) J : 



q=l 



n-1 

^(-1)"+PL> 
p=l 

(0, - C) 



din-ih 

di 



dipk + D 



A; 



«l«n-l 



n-1 



d. 



-''«n-l«l 

diph 



T'Ti — l in — 1 


^n — 1 




djiq 




dii 


din-lin~l 


din- 


dipi„-i 


dip 



+ D 



3 
k 



iCk - c; ) 



D. 
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It is easily seen that in the last expression, only the last term remains. This is the right-hand side 
of (EH. □ 

Corollary 5.2. The vector u^i-^^... defined for distinct positive integers ii, ■■ ■ ,in-i,in by 

d-^- j5D with ^4' 'i'S 0, unit vector, i.e. 

Proof. Using Lemma l5. II in the special case j = k{= in), we have 
II 11^ — ^ 



n-1 /■_/■* 2 



n 

1=1 





diiii 


diii„_i 












din-i'h 


din-lin-1 




diiii 


dilin 






dinh 


dinin 





n-1 

n 

1=1 



A* 





diiii 


diiin~\ 




diiii 


dilin 


i 


din-lh 


di„-ii„-i 




dinh 


din^n 



diiii 



-^*n-l«l 



diii„-i 



^n — 1 ^n — 1 



□ 



We are now able to apply the collision laws defined by Theorem 13.11 to the two-soliton colli- 
sion in which soliton Q{i-^^... ^i^_-i^j}^k overtakes soliton [cf. the definition H4-16() ]. Here, 
ii, - ■ ■ ,in-i,j, k are distinct positive integers. 



Proposition 5.3. The two-soliton collision in which overtakes Qji^,... changes 

these solitons to q{i^^... ^i^_-^]^k '^'^d ^{i^^... ■i„_-^^k},j! c-^ shown in Fig.\^ According to Theorem XS.li 
this is equivalent to the following set of equalities: 



Cj - Ck 



- Q 



1 + 



10 -CP 



■"{n,-,in-i},i ■ '"'{ii,- ,in-i},k 



(5-2a) 
(5-2b) 



/-* 

. , - p'/'{h,-,in-l,i),k-<l>{h,-,i„-i},k2i I^.<ur. . 1, 
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- Cj / t 



(5-3) 



{il,--- ,i„-l,k),j <P{n,--- ,in-l}.J 



/-^ 

^ ^ 

Ck - q 



u 



{«!,■■■ ,in-l},j 



Ck - Q 

Ck - Ci 



Us 



(5-4) 



Proof. Throughout this proof, we employ the following notation in order to express determinants 
compactly: 



J1,J2,--- ,Jl 



djiki 'djik'z 
dj2ki dj2k2 



djiki 



d 



djih 

dj2ki 



To prove (|5-2j) . we first rewrite the left-hand side of (|5-2a|) as 



Cj - Ck 



- ct 



d 



-1,3 



*i ) ■ ■ ■ ) 1 ) J 



) ^n— 1 ) ^ 
^1 ) ■ ■ ■ ) ^n— 1 ) ^ 



(5-5) 



using the definition of 4'{i'i.— l|4-14() . Obviously, the right-hand side of (jS-Sf) is symmetric 

with respect to interchange of j and k. It follows that the equality ()5-2ap holds. 

Next, we prove the equality (|5-2b|) . Using the definition oiui^i^ ... ()4-15|) . and Lemma l^TTl 

we obtain 



{il,--- ,i„_l>,i+'f'{n,--- ■^n-l},'! 



Cfc - C,* Vr (C* - C*j){Cii - Ck) 

^\ {Cn-qm,-Ck) 



I'll ' ' ' 1 ^n— li J 



• (5-6) 



^li ■ ' ' 5 1 

Multiplying (|5-6() by its complex conjugate on each side, we obtain, with the help of H4-14() 



"^{^ir" i^n— "^-{ii,--- , in— ill's 



a - c: 



d 



n, • • 
n, • • 



(o-c;)(a-a) ^mi,-- 



; ^n— 1 ) ^ 



^1 ; ■ ■ ■ ; ^n— 1 ) J 



^1 ) ' ■ ■ 1 ^n— 1 1 ] 



A; 



^1 ) ■ ■ ■ ) ^n — 1 ; 



Then, we can rewrite the right-hand side of H5-2b|) as 



- Ck 

- a 



[ (o-c;)(a-a) 
I 10 -ap 



{il,--- ,i„_i},i ■ "{i^, 
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Cj - Ck 



^1 ; ■ ■ ■ 5 ^n— 1 ; J 



^1 ) ' ■ ■ 7 ^n— 1 5 ^ 



*1 ; ■ ■ ■ 5 1 ; J 



«1, 



-1,J 



^1 ) ■ ■ ■ 7 ^n— 1 ) ^ 
^1 ) ■ ■ ■ 1 ^n— 1 5 ^ 



(5-7) 



Here, thanks to the Jacobi formula for determinants, we have 





1 ^n— 1 




ii, ■ ■ 


) 1 


i 


n, • • 


) ^n— 1 


) 


h, ■ ■ 


) ^n— 1 





■ 5 in~li k 

■ 5 ^71— 1 ) k 

I ^n— 1 ) J) ^ 



• • • 



-1, 



n, 



) ^n— 1 



(5-8) 

We remark that, 



Thus, 1)5 -7^ is equal to (jS-Sp . This completes the proof of the equality H5-2b|) 

To prove the equality (|5-3jl . we need to extend the Jacobi formula 
although the matrix elements du here are given by H4-8|) . the two sides of ()5-8p are equal as a 
polynomial for general elements du. Therefore, maintaining the validity of (|5-8() . we can replace 
the columns with index k by columns consisting of the vectors , 

di 



,Ui,,_i,Uj or life: 



^1, 



^1,3 



diiii 



d.^ 



Hill 



-in 



di^ 



-l«n-l 



dki„ 



Ui 



Uk 



d 



I 1 5 



Ml 



^1 ) ' ' ' ) ^n— 1 
^1 ! ' ' ' ) ^n— 1 



djii 
dkii 



^n — l^n — l 

dki„-i 



dkj 



u 



«n-l 



(5-9) 



We rewrite the right-hand side of ()5-3() using H4-14() . H4-15() and H5-6() (with j ^ k) as 

/-^ A* 



_ Q'Y^{h,- ,in-i},k _ ^* \^'"{nr--,in-i},fc ■ ■"■{^l,...,^„_l},j 



X 11 



{«!,■■■ ,i„~i},j 



/-* A* n— 1 A* 

{n.-,i„-i,i},fc X — Yl — — X 



Cfc di Cfc 



X < 



^1 5 ■ ■ ■ ; ^n— 1 ) J 



dkil 







• • • , 


«n-l, j 


diiin-i 




1. 

^n — 1 ^n — 1 


Uj 


dki„-i 


Uk 



^1 5 ' ■ ■ ; ^n— 1 
^1 ) ' ■ ■ ; ^n— 1 



^1 ) ■ ■ ■ ) in—li ^ 
^1 ) ■ ■ ■ ) ^n— 1 ) J 
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-111 



diii 



din-lin-l 



Ui 



u 



d 



Owing to the extended Jacobi formula (|5-9|) . this is equal to u^^-^ ... ^i^_j^jy^k [cf- the definition (|4-15|) ]. 
Now, the proof of the equality (|5v^|) is complete. The proof of the equality (|5-4|) is obtained by 
interchanging j and k in the proof of H5-3() . □ 



Proposition 15.31 is applicable to two-soliton collisions satisfying the following condition: 

• Before the collision, the set of subscripts in the left-hand soliton's { } is equal to the entire 
set of subscripts of the right-hand soliton: {ii,--- ,in-i, j} = {h, ■ " )^n~i}U{j}. 

We also note the following properties (cf. Fig. ISJ: 

• After the collision, the set of subscripts in the left-hand soliton's { } is still equal to the entire 
set of subscripts of the right-hand soliton: {ii, • • • , in-i, k} = {ii, ■ ■ ■ , in-i} U {A;}. 

• The entire set of subscripts of the left-hand soliton is unchanged: {ii, • • • j} U {k} = 
{ii,--- ,in-i,k}U {j}. 

• The set of subscripts in the right-hand soliton's { } is unchanged: {ii, • • • , in-i} = {^i, • ' ' j in-i}- 

• The overtaken soliton's number is removed from the overtaking soliton's { }, while the over- 
taking soliton's number is added to the overtaken soliton's { }. 



Q{il,...,in-i,k},j 




; interaction 
range 




Fig. 3. Two-soliton collision in the presence of other solitons. 
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We are now able to state the main result of this paper. 

Theorem 5.4. An N-soliton collision in the Manakov model can he factorized into a nonlinear 
superposition of (^) pair collisions in arbitrary order. 

Proof. According to the asymptotic behavior of the A^-soliton solution as t ^ — oo (see Fig. I^J, 
solitons-A^, • • • , 1 are initially distributed along the 

- ,Af-l},Af) ■ ■ ■ 1 ,n-l},ni " " " ' 9{ },1- 

We take this initial state as the point of departure and assume that the solitons collide pairwise in 
a given order. Then, a pair collision takes place (^) = ^^^^^^ times. What will the final state be 
under this assumption? To answer this, we note the following two points: 

• The set of subscripts in each soliton's { } is always equal to the entire set of subscripts of the 
next soliton to the right. This ensures that Proposition ESI is applicable to every pair collision. 

• Soliton-n will overtake solitons-1, • • • , n — 1 and will be overtaken by solitons-n + 1, ■ ■ ■ , N. 
We see that, regardless of the order of the pair collisions, solitons-1, ■ ■ ■ , N are finally distributed 
along the x-axis as 

Q{2,--- ,N},1: ■ ■ ■ ) Q{n+1,--- ,N},n^ ' " " > 9{ },7V- 

This final state is exactly the same as the asymptotic behavior of the A^-soliton solution in the 
t — ^ +00 hmit (see Fig. 

Q.E.D. 

Remark 1. Theorem 15.41 together with (|H-l()j) . demonstrates that the initial state of solitons 
uniquely determines its final state. That is, an A^-soliton collision described by Proposition 14.21 
defines a mapping from the initial state to the final state. This fact is not evident from Proposi- 
tion O 

Remark 2. We can now prove that the mapping p-lUj) satisfies the property (b) or, equivalently 
the Yang-Baxter property (b'), as stated in the introduction. It is sufficient to verify the following: 

The initial polarization vectors of N solitons (wjx^... ^iv-i},Af; ' ' ' 7'^{},i) given in 
Proposition 14.21 can be made to coincide with any combination of N unit vectors in 
C" by appropriately varying the bare polarization vectors (itAr, • • • 
This is easily verified if we consider the order of the pair collisions such that every soliton ex- 
periences the bare polarization U{},j(= Uj) once and inductively use the fact that the mapping 
U2 I— > i*{i},2 ill Theorem 13. II is a bijection on the unit sphere. 

Remark 3. The validity of the Yang-Baxter property allows us to extract a new "set-theoretical" 
solution to the quantum Yang-Baxter equation (cf. Refs. 15) and 1^). To be more specific, re- 
garding (|3-10() as the mapping (u{i} 2,1^1) 1— > (w2iW{2}^i), we obtain a nontrivial solution to the 
parameter-dependent Yang-Baxter equation for mappings that act on the direct product of two 
(complex) unit vectors. Naturally, we can extend this solution further by adding information con- 
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cerning the center positions of solitons. However, this extension is not very intriguing, because the 
net change wrought by the mapping does not depend on the added information. 

§6. Concluding remarks 

In this paper, we have investigated sohton colhsions in the Manakov model for the general case 
of m components using a straightforward approach. We first derived the general A^-soliton 

solution of the Manakov model from that of the matrix NLS equation ()2-2|) through a simple 
reduction.® We considered the limits t ^oo for the N = 2 case and obtained the collision laws 
of two solitons in the Manakov model. Next, we considered the same limits for the case of general 
N and obtained the asymptotic behavior of the A^-soliton solution. We were able to diagram the 
asymptotic behavior in the simplest way in terms of the quantity defined by H4-16|) 

(see Fig. Taking advantage of this, we proved with a simple combinatorial treatment that an 
A^-soliton collision in the Manakov model can be factorized into a nonlinear superposition of (■^) 
pair collisions in arbitrary order. This clears up the longtime misunderstanding that multi-particle 
effects exist in the Manakov model. 

This result is far from trivial in the m > 2 case. In the m = 1 case (scalar NLS), all the 
soliton parameters that play an essential role in the collision laws (in the notation of this paper, 
Cij C2; ■ ■ ■ ) Cn) are invariant in time. A pair collision results only in a displacement of the soliton 
centers and a shift of the phases, which will not change the effects of future pair collisions. Thus, a 
superposition of (^) pair collisions gives the same results for every order of the pair collisions. It 
is not difficult to prove in this case that an A^-soliton collision reduces to a pair collision .^'^'C^ 
In contrast, in the m > 2 case, a pair collision results in a change of the polarization vectors. This 
changes the effects of future pair collisions completely. Therefore, it was not obvious before the 
present work that a nonlinear superposition of (^) pair collisions gives the same results for every 
order of the pair collisions or that it exactly coincides with an A^-soliton collision. The key to 
proving these facts is a highly nontrivial relation among determinants and extended determinants, 
given in Lemma [5. 11 This implies the possibility that some new relations similar to Lemma IS.ll can 
be obtained through the investigation of soliton collisions in multi-component integrable systems. 

Very recently, Steiglitz and coworkerj^'E3 proposed that soliton collisions in the Manakov 
model can be utilized to carry out any computation with beams in a nonlinear optical medium 
(see Ref. I38|) for the experimental foundations). We believe that the results obtained in this paper 
will be useful to refine and reinforce their interesting idea. In particular. Theorem 15.41 supports*^ 
their hypothesis on the pairwise nature of soliton collisions and Proposition 14.21 together with 
Proposition 15.31 provides a hint on how to design optical logic operations more simply and reliably 
than the method proposed in Ref. 37). 



*' To be precise, we need to extend our result to the case in which some solitons propagate at the same velocity 
and do not collide with one another. This is beyond the scope of this paper and is left to the reader as a future 
problem. 
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